Emission line studies of young radio sources 
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Abstract 

£^ . We demonstrate the efficiency of high quality optical spectroscopic observations 

of two compact radio sources, PKS 1549-79 and PKS 1345+12, as a probe of the 
■ kinematics and physical conditions in the circumnuclear gas in the early stages of 

radio source evolution. We outline a schematic model for PKS 1345+12 based on 

£T) \ the model for PKS 1549-79 proposed by Tadhunter et al. (2001) in which the young 

radio source is expanding through a dense and dusty enshrouding cocoon, sweeping 

^ ■ material out of the circumnuclear regions. 
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C$ 1 1 Introduction 



Gigahertz-Peaked Spectrum Radio Sources (GPS: D < 1 kpc) and the larger 
Compact Steep Spectrum Radio Sources (CSS: D < 15 kpc) account for a 
significant fraction of the radio source population (~ 40%) though their nature 
is not fully understood (see O'Dea 1998 and references therein). Currently, we 
believe they are young radio sources (Fanti et al. 1995) supported by estimates 
of dynamical ages: td yn ~ 10 2 - 10 3 years (Owsianik et al. 1998); and radio 
spectral ages: t sp < 10 4 years (Murgia et al. 1999). This is in preference to the 
frustration scenario where the ISM is so dense, the radio jets cannot escape 
and the radio source remains confined and frustrated for its entire lifetime 
(van Breugel 1984). 
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Fig. 1. Schematic model for PKS 1549-79 proposed by Tadhunter et al. (2001). 

If compact radio sources are young radio sources, we will observe them rela- 
tively recently after the event (s) which triggered the activity (e.g. a merger; 
Heckman et al. 1986). Indeed, many compact radio sources exhibit features 
attributed to mergers such as double nuclei, tidal tails, arcs of emission and 
distorted isophotes (e.g. Stanghellini et al. 1993). Hence, the host galaxy will 
still retain the dense ISM injected by the merger. The young radio jets will 
also be relatively small and so will readily interact with the dense ISM. 

In these proceedings we demonstrate the efficiency of high quality optical 
spectroscopic observations as a probe of the kinematics and physical conditions 
of the ISM in the host galaxies of radio sources, particularly in the early 
stages of radio source evolution. We discuss observations of two compact radio 
sources, the flat spectrum radio source PKS 1549-79 and the GPS source 
PKS 1345+12, and outline a schematic model for PKS 1345+12. For detailed 
discussions see Tadhunter et al. (2001) and Holt et al. (2002, 2003). 



2 PKS 1549-79 

From low resolution spectra of the compact flat-spectrum radio source PKS 
1549-79, Tadhunter et al. (2001) report: 

• highly broadened forbidden emission lines (FWHM ~ 1350 km s _1 ). 

• 2 distinct redshift systems: the high ionisation forbidden lines (e.g. [O III]) 
are blueshifted by ~ 600 km s _1 with respect to the low ionisation lines 
(e.g. [O II]) and HI 21 cm absorption. 
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Fig. 2. Left: Radial velocity profile of the extended gaseous halo. Key: *, O and 
□ : highly extended [O II] emission (narrowest component); •: HI 21 cm absorption 
(Mirabel 1989); A: 3 components of [O III] in the nuclues, from right). 

Tadhunter et al. (2001) interpret the blueshifted emission lines as material in 
outflow. Due to the small scale radio source (~ 300 pc), the radio morphology 
(core-jet) and the flat radio spectrum, they proposed a schematic model in 
which the young radio source is enshrouded by a cocoon of dense gas and 
dust completely obscuring the central QSO nucleus. The young radio jets are 
expanding and sweeping material out of the nuclear regions with the direction 
of propagation close to the observer's line-of-sight (see Fig 1). 



3 PKS 1345+12 



We present here a summary of results from intermediate resolution (~ 4A) 
spectra over a wide spectral range (~ 4500A) from ISIS on the 4.2m WHT on 
La Palma. For a detailed discussion see Holt et al. (2003). 

From modelling of the highly extended [O II]AA3726,3729 emission line (~ 
20 kpc) using between 1 and 3 Gaussian components, we believe that the 
narrowest of these represents the galaxy rest frame: 

• It is the only component observed across the entire spatial range; 

• It is consistent with the deep HI absorption (Mirabel 1989; Morganti et al. 
2002) (see Fig 2) and the stellar absorption lines (Grandi 1977); 

• It is the least kinematically disturbed component (FWHM ~ 340 km s _1 ); 
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• It has a small velocity amplitude (A ve i < 250 km s _1 ; see Fig 2) consistent 
with gravitational motions (Tadhunter et al. 1989); 

• It is the only component consistently observed in all emission lines, irre- 
spective of the model required to reproduce the emission line profile. 

In the nuclear aperture, the emission lines are broad with strong blue asym- 
metries. It is essential to model all emission lines with 3 Gaussian components 
(see Fig 2, also Holt et al. 2003). The narrowest of these is consistent with the 
narrow component of [O II] and the rest frame of the galaxy. The intermediate 
and broad components are blueshifted by up to ~ 2000 km s _1 with respect to 
this narrow component. Due to large reddening in the nucleus (see below), we 
interpret this as material in outflow. However, one model does not reproduce 
all emission lines - different velocity widths and shifts for the broad and inter- 
mediate components are required for [O I]AA6300,6363 and [S II]AA6716,6731 
(see Holt et al. 2003). 

We estimated the density in the nucleus using the density diagnostic, the [S 
II]AA6716,6731 doublet. The intermediate and broad components have densi- 
ties consistent with the high density limit. By varying the fit parameters, we 
obtained lower limits on on the density of n e > 5300 cm -3 and n e > 4200 
cm~ 3 for the intermediate and broad components respectively. The narrow 
component is consistent with the low density limit (n e < 150 cm -3 ). 

We have investigated reddening in the nucleus of PKS 1345+12 using three 
independent techniques, each assuming a simple foreground screen for inter- 
stellar extinction and the Seaton (1979) extinction law. 

• B aimer line ratios. The Hct/[N II] blend is highly complex. We measure 
Ha/H/3 ratios of 3.32 ± 0.33, 5.25 ± 0.28 and 18.81 ± 4.74 corresponding 
to E(B-V) values of 0.06 ± 0.05, 0.42 ± 0.10 and 1.44 ± 0.50 for the narrow, 
intermediate and broad components respectively. By varying the model fit- 
ting parameters, we estimate a lower limit on the reddening in the broadest 
component of E(B-V) > 0.92. 

• Comparison with the infra-red. Veilleux et al. (1997) present measurements 
of Pact. The broadest component is consistent in velocity width and shift 
with the broadest component in [O III]. By varying the quoted flux of Paa 
by a factor of 2, to account for possible slit loss differences between the 
observations, Paa/H/3 gives a range of E(B-V) values: 1.60 < E(B-V) < 2.00. 
This is consistent, within the errors, with the measured value for Ha/H/5. 

• The nebular continuum. Before modelling the faint lines in the nucleus, we 
subtracted both the nebular (see Dickson et al. 1995) and stellar contin- 
uum. If zero reddening is assumed for all components, a large discontinuity, 
corresponding to the Balmer edge at ~ 3645 A, remains when subtracting 
the nebular continuum from the spectrum (see Fig 7 in Holt et al. 2003). 
However, if the measured E(B-V) values are assumed, the discontinuity dis- 
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Fig. 3. Schematic model for PKS 1345+12. 

appears. The strength of this discontinuity is particularly sensitive to the 
amount of flux in the intermediate component and we estimate a lower limit 
of E(B-V) > 0.3 for the intermediate component. 

The nucleus of PKS 1345+12 is highly reddened, with reddening increasing 
with broadness of component. 

We estimate an upper limit on the mass of line emitting gas of < 10 6 M 
(from the density and reddening-corrected FL3 luminosity) which is consistent 
with the radio source being young. 



4 The model 



Combining these results, we have developed the model from Tadhunter et al. 
(2001). Fig 3 shows the schematic model which retains the small scale ra- 
dio jets expanding out through an enshrouding cocoon giving rise to emission 
line outflows. New features include stratification of the ISM - all emission 
lines in the nucleus require 3 Gaussian components though the velocity width 
and shift of the intermediate and broad components varies between emission 
lines. However, gradients (e.g. density, velocity, ionisation potential) must ex- 
ist across the regions emitting the intermediate and broad components. The 
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narrow component remains consistent in all emission lines and represents the 
relatively quiescent halo. 'Strata' positioning is determined by the reddening 
- the reddening increases with broadness of component and so the broadest 
component must originate from the region closest to the nucleus. 
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